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A dual foliation treatment of General Relativity is presented. The basic idea of the construction 
is to consider two foliations of a spacetime by spacelike hypersurfaces and relate the two geometries. 
The treatment is expected to be useful in various situations, and in particular whenever one would 
like to compare objects represented in different coordinates. Potential examples include the con¬ 
struction of initial data and the study of trapped tubes. It is common for studies in mathematical 
relativity to employ a double-null gauge. In such studies local well-posedness is treated by referring 
back, for example, to the generalized harmonic formulation, global properties of solutions being 
treated in a separate formalism. As a hrst application of the dual foliation formulation we find 
that one can in fact obtain local well-posedness in the double-null coordinates directly, which should 
allow their use in numerical relativity with standard methods. With due care it is expected that 
practically any coordinates can be used with this approach. 


I. INTRODUCTION 

For their consideration as an initial value problem the 
field equations of General Relativity are typically split 
into a set of evolution and constraint equations. This is 
done by introducing coordinates The level 

sets of the time coordinate t are taken to be spacelike hy¬ 
persurfaces which foliate the spacetime. The unit normal 
to the hypersurfaces is used to 3 -|- 1 decompose the field 
equations in the natural way. This results in the vacuum 
field equations in the textbook form 

^tUij — ‘ZcyKij Gjn'yij , 

dtKij = —DiDjU + a,\Rij — 2K’^iKjk + KKij] + CpKij , 

H = R- +K^ = 0, 

M, = D,K = 0. (1) 

Given two sets of observers, one associated with , an¬ 
other with coordinates XT- = {T,X-), what is the re¬ 
lationship between the solutions as expressed in each 
in the 3-1-1 picture? Unfortunately a clear presenta¬ 
tion of the resulting formalism is not readily available, 
despite being straightforwardly obtained by space-time 
decomposition of the four-dimensional Jacobians JT^ = 
dXT/dx^. The first aim of this work is to give just such a 
presentation, which can be found in section im This dual 
foliation approach will be useful in numerical relativity, 
where one expects it will help in the construction of ini¬ 
tial data and in the comparison of solutions constructed 
with different choices of lapse a and shift /3*. 

Amongst the most powerful machinery in mathemati¬ 
cal relativity is that of the double-null coordinates. With 
this choice the field equations exhibit a particular struc¬ 
ture that allows the demonstration of the stability of 
Minkowski spacetime Q , and that a certain special class 
of vacuum initial data will collapse to form blackholes . 
One would thus like to use these coordinates in numerical 
relativity, preferably with standard methods. A number 
of applications present themselves; the conjectured insta¬ 
bility of Gauchy horizons the propagation of weak- 


null siMularities Q , and the critical formation of black- 
holes [8|. Unfortunately from this perspective the proofs 
of these impressive results employ a different formalism 
for long-term results and local existence. This is a serious 
bugbear because, as painful experience has taught, a nec¬ 
essary condition for any numerical method to converge is 
that the PDF problem is locally well-posed. Therefore 
the second aim of this work is to find such a formulation. 
In section Hill this is attempted in a direct way. The nor¬ 
mal approach is to modify the equations by introducing 
new constraints, coupled in a particular way to the gauge 
choice, and insodoing uncover, say, a strongly hyperbolic 
formulation. But we find in a pure gauge analysis that 
the standard form of the double-null coordinates are orfiy 
weakly hyperbolic, and so can not be used in this way [9|- 
[TTj| . With appropriate alterations, there may be such a 
straightforward formulation, but because a preferred di¬ 
rection is singled out the construction will in any case be 
complicated. We thus abandon the search. 

In section IIVI we summarize the first order reduc¬ 
tion [l^ of t he g eneralized harmonic gauge (GHG) for¬ 
mulation |13M13 employed in the numerical relativity 
codes SpEG [ig and bamps [IMl- We use the dual 
foliation formalism to derive equations of motion for the 
Jacobian that maps from generalized harmonic to double¬ 
null coordinates. These equations are minimally coupled 
to the field equations, and so their hyperbolicity may be 
treated easily. Indeed the Jacobians satisfy a set of non¬ 
linear advection equations, and so are hyperbolic. We 
may consider the full set of fields to be solved for as the 
GHG system with the Jacobians tacked on. We can sub¬ 
sequently change independent variables from x^ -T XT, 
with XT the double-null coordinates. The punchline is 
that since the system has at most first derivatives, we can 
do so without generating any derivatives of the evolved 
Jacobians. Therefore the PDE properties of the system 
are unaffected and we end up with a formulation that 
is symmetric hyperbolic directly in the double-null coor¬ 
dinates. Weaker notions of hyperbolicity are also con¬ 
sidered. Concluding remarks are collected in section IVl 
Geometric units are used throughout. 
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II. THE DUAL FOLIATION FORMALISM 

In this section we work in the intersection of two co¬ 
ordinate patches = (t,x^) and AF = (T, A-), related 
by the Jacobian = dXi^/dx^. The two time coor¬ 
dinates define distinct foliations of the spacetime, and 
with them different notions of spacelike tensors, intrin¬ 
sic and extrinsic curvatures. These quantities are related 
in the natural way with a 3 -I- 1 split of the Jacobian. 
Consequently the form of the gravitational field equa¬ 
tions in each foliation is related. Throughout latin in¬ 
dices a, b, c, d, e will be abstract. Greek indices stand for 
those in coordinates x^, or if underlined in AF. Simi¬ 
larly latin indices i, j^k,l,m,p stand for spatial compo¬ 
nents in a;^, and when underlined for spatial components 
in Alt. Indices n and v denote contraction in that slot 
with the vectors n“ and v°‘ respectively. The summation 
convention is always employed. 


A. Coordinate freedom 

Coordinate change under a 3 -I-1 decomposition: Con¬ 
sider two sets of coordinates x^ and AF defined in the 
same region of spacetime. Each of the two time coordi¬ 
nates t and T naturally defines a foliation of the space- 
time which we will refer to as the lower case and upper 
case foliations respectively. In the lower case foliation 
we define the standard lapse, normal vector, time vector, 
projection operator, and shift vector, 

a = {-Vat , 

CVat = 1, T“b = 5\ + n^Ub , 

I3a=±\tb, P^ = -an‘^VaX\ ( 2 ) 

With both indices downstairs the projection operator 
becomes the natural induced metric 7 o& on the lower case 
foliation. The covariant derivative associated with ^ab is 
denoted by D with connection T. The same definitions 
are made in the upper case foliation, 

A= (-VaTV“r)-^ , N‘^ = -AV‘^T, 

T'^VaT = 1, -t A“ Ab , 

Ba = , B^- = -AN'^VaX^-. (3) 

The covariant derivative associated with is denoted 
by with connection 

The Lorentz factor and boost vector: The unit normal 
vectors of the upper and lower case foliations are related 

by 


A“ = IT(n“-kz;“), (4) 

where we have defined the Lorentz factor W and lower 
case boost vector Va, 


= ^l-\Nb. (5) 


TABLE I: A summary of the definitions of the various met¬ 
rics, time reduction variables the relationship between them. 
The fourth column gives the object used as a time reduction 
variable when employing the given metric. The final column 
gives states the equation numbers relating the curvature of 
the given metric with that of the others. ‘GCM’ stands for 
the Gauss-Codazzi-Mainardi equations. 


Metric 

Connection 

Defn. 

Time der. 

Curv. 

Note 

9ab 

v,(4)r 



^‘^'^Rab 


'yab 

D,r 

-L 9ah 

Kab 

^ab 

‘GCM’ 


D,(E 

T 

IKa6 

Rab 

(|33l) 

Jab 

(SI)JJ (N)p 


^^^Kab 


(|37ll 

Since 

the normal 

vectors 

have unit 

magnitude the 


Lorentz factor and boost vector satisfy, 

W = = , IE > 1 > Y^ViVj = v'^ . (6) 

vl — ViV'^ 

This is simply the requirement that the upper and lower 
case normal observers travel subluminally relative to one 
another. Observe that we also have, 


= IE(A“ -b E“), (7) 


with the upper case boost vector, 

Ea = ^^^^l-^arib = {W-^ - W)na - WVa , (8) 


so that there is a natural reciprocity between the coordi¬ 
nate systems as expected. We thus also obtain. 


W = 


1 

v'l - W-' 


( 9 ) 


Provided the spatial boost vector, a vector Sa which is 
spacelike in the upper case foliation, SaX^- = 0, can be 
reconstructed from its projection into the lower case fo¬ 
liation with, 

+ ( 10 ) 


and obviously a similar result holds for all tensor va¬ 
lences. Therefore we may restrict our attention of upper 
case spacelike tensors to their projections into the lower 
case foliation, and thus look only at the spatial compo¬ 
nents in the lower case tensor basis. 

Decomposition of the Jacobian ; Upon 3-1-1 de¬ 
composition the Jacobian = dXiJjdx^ can be writ¬ 
ten, 


J-a A„ = -W , n“ = TT^ , 

J^,N^=Wv,, Jh = 4>^^. ( 11 ) 

The components tt- are given in terms of the upper case 
lapse, shift and boost vectors by, 

= WV^-- WA-^B^-, 


W = -{N'^na) 


V, 


( 12 ) 
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although we rarely find that this is the most convenient 
form for the expression. In matrix form we therefore have 
the representation, 

j^( A-^W{a-P'^Vi) a7r2- +\ 

V -A-^Wv, ^ ’ 

Note that by introducing the variables A~^Wvi and (j^i 
to replace first order spatial derivatives of the coordinates 
we have effectively introduced reduction constraints, 

D[,{A-^Wv,}) = 0 , = 0 , ( 14 ) 

which we will refer to as the hypersurface (orthogonality) 
constraints. Here and in what follows one must be care¬ 
ful to note that the upper case underlined index is to be 
treated as a simple label rather than an open slot when 
working in the lower case coordinates. It is straightfor¬ 
wardly checked that the transformation, 

a-^W{A- 

-a-^WV, ' y ’ 

with the various auxiliary quantities defined in the nat¬ 
ural way, is indeed the inverse transformation. 

Time development of the Jaeobian: By the equality 
of mixed partials we have the Hamilton-Jacobi like equa¬ 
tions, 

dt{A-^Wv,) = -D,[a{A-^W)\ + Lp[A-^Wvi), 

= + (16) 

for the components J—i and J\. These expressions hold 
regardless of the upper case coordinate choice, but the re¬ 
maining four components can be determined only once a 
particular coordinate choice is known. Perhaps the sim¬ 
plest useful example is the generalized harmonic gauge 
choice DX— = H—, which results in, 

dt{A-^W) = a{A-^W){K + E) - D\aA-^Wv^) 

+ Cp{A-^W), 

dtTT^ = {a<lAj) + aEi- + , ( 17 ) 

where the gauge source functions are decomposed as E = 
{A/W) E[— and E- = -El-. The lower case extrinsic 
curvature is defined by, 

Kab = — T'^aVcU-b , (18) 

and likewise in the upper case foliation, except that as 
elsewhere we append a label N. On a given spacetime 
with coordinates the system (UHl), (HZl) can be viewed 
as a simple first order reduction of the four wave equa¬ 
tions 'C\Xlt- = iJii. Other choices will be more conve¬ 
niently expressed once the relationship between the two 
induced geometries are known. 


Equations of motion for projected upper case objects: 
The equations of motion of the projection of upper case 
spacelike tensors Sa and Sab projected into the lower case 
foliation are, 

Q; 

Hn Si ^(^l3—a.v)Si ; 

O' 

dtSij = Sij -\- i^(^^ — oiv')Sij , ( 19 ) 

for vectors and symmetric tensors respectively. Similar 
expressions hold for arbitrary valences but will not be 
used in what follows. The projected upper case induced 
metric defined by ga& = is, 

gy = ''"Vy = 7y + W^ViVj . (20) 

Note that as the sum of a symmetric positive definite 
matrix and and semi-positive definite combination of the 
boost vector, the projected upper case metric is itself 
positive definite and thus invertible, and can be consid¬ 
ered a metric on leaves of the lower case foliation in its 
own right, if we so wish. When doing so we will refer to 
this object as the boost metric, and denote the covariant 
derivative by D with connection €. The boost metric has 
inverse, 

(g-i)y = yj _ ^ (21) 

by the Sherman-Morrison formula. We now aim to relate 
the geometry of the upper and lower case foliations in 
terms of the lower case normal, Lorentz factor and boost 
vector. A summary of the relationships between the four 
different metrics gab, Jab, ^^^7ab and gab is given in Table ID 
Connections and curvatures: The Christoffel symbol 
associated with the upper case induced metric is given 
by the standard expression, 

, (22) 

which holds in arbitrary coordinates, and where here and 
in what follows we use the projection operator without 
indices to denote the projection on every open slot. The 
Christoffel symbol associated with the lower case induced 
metric is defined similarly. By the argument around 
equation (flUl) we need only compute the projection of 
the upper case Christoffel into the lower case foliation. 
Using (EiH) and (1^ we Hnd that, 

g"*, gp, - 2 w2 g"' 

+ g™, gP, Kmp - 2 

+ W^g^^ ai + 2W^a-^ g^ 

-I- W^a~^ g^i ViVj dn/3'' - W^v^ViVj Ina , 

(23) 

where here we use an index n to denote contraction with 
the lower case unit normal vector , and where the accel¬ 
eration of lower case Eulerian observers is Oi = Dihia. 
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The upper case induced Ricci tensor can be computed 
from, 


+ , ( 24 ) 

and likewise for the lower case curvature. The rela¬ 
tionships between the upper and lower case connections 
above (1^^ can be used to compute the relationship be¬ 
tween the two Ricci curvatures by brute force, but it is 
more convenient to use the Gauss-Codazzi equations, as 
described momentarily. The upper case extrinsic curva¬ 
ture projected into the lower case foliation is, 

= WK,, - - WAf^iVj) = WiK,j - , 

( 25 ) 

where here we also define IKy which, for lack of a better 
name, we call the boost extrinsic curvature. Note that 
we define the trace IK = (g~^)®'lKy in a nonstandard way 
by using the inverse boost metric. The projected upper 
case acceleration is, 

lki = Ai= DiilnA) + DjihiA) + £„(lnA)^ . 

( 26 ) 

It is most convenient to express the various equations 
in terms of Ai rather than using this expression, as we 
wish to suppress the explicit appearance of a £„(lnA) 
contribution. 

Constraints under the coordinate change: Comparing 
the Hamiltonian and momentum constraints in each fo¬ 
liation, we find that, 

= W^H - 2W^My , 

= WM, + 2W^M^Vi - W^Hvi. ( 27 ) 

An index v denotes contraction with the boost vector ti“. 

It is thus obvious that the full set of constraints will 
be satisfied in the lower case foliation if and only if 
they are satisfied in the upper case foliation. Expand¬ 
ing out the upper and lower case constraints and using 
projected upper case extrinsic curvature (EH) in combi¬ 
nation with dini), we easily find the relationship between 
the two spatial Ricci scalars. As it stands, equation (l?7l) 
is really of a purely geometric nature and independent 
of the gravitational field equations, where we think of 
the symbols as mere shorthands according to CD, or the 
upper case foliation equivalent. 

Electric and Magnetic decomposition of the Weyl ten¬ 
sor: Especially in General Relativity the decomposition 
of the Weyl tensor Wabcd into two spatial tracefree ten¬ 
sors, Electric and Magnetic parts, 

Eab = n^n‘^WacM , Bab = n^n‘^ *Wacbd , ( 28 ) 

has special significance in encoding the propagating de¬ 
grees of freedom of the gravitational field. The dual Weyl 


tensor here is * Wabcd = cdWgfab- Evidently this de¬ 
composition is foliation dependent, because of the pres¬ 
ence of the normal vector n“. The relationship between 
the two decompositions is given by, 

= {2W^ - l)Ay - 

<'^’Ry = W^Hy - lE^e'y , (29) 

where Cbcd = n°^£abcd stands for the lower case spatial vol¬ 
ume form. Furthermore this equation shows trivially that 
changes of coordinates can not create nor destroy gravita¬ 
tional waves. Since in vacuum the electric and magnetic 
parts also satisfy a closed evolution system , it is also 
clear that if we are given initial data with vanishing Ay 
and Bij this will be the case once and for all. The re¬ 
lationship (1291) holds in general, but using the vacuum 
Einstein equations, we have that, 

Eab = - (K^Kbcf^ + KKJ! , (30) 

and likewise for the upper case Electric part. We can 
thus relate the tracefree part of the upper case and lower 
case spatial Ricci tensors as we did for the Ricci scalars, 
namely by expanding out the projected upper case ex¬ 
trinsic curvature with (ESI) and using (fTUl) to obtain the 
non-spatial components. 

Boost metric connection and curvature: We would 
like to have the equations of motion for upper case ob¬ 
jects in the lower case coordinates. But as it is more 
convenient to work with spatial tensors in the lower case 
coordinates we instead work with the boost metric and 
boost extrinsic curvature (gy, Ky) to obtain the desired 
results. The time derivative of the boost vector is conve¬ 
niently encoded as, 

dt(Wvi) = a W[A, - A„u, - D, ln(a W)] + Cp{Wv ^). 

( 31 ) 


The relationship between the connection of the boost 
metric and the lower case spatial curvature is C^y = 
r'^y , with 




2 u Ey'yij 


(32) 


This result can be obtained from the standard expression, 
see Ch. 7 of [2l|, since the lower case spatial metric and 
boost metric act in the same tangent space; notice that 
this is not the case when we try to relate the upper and 
lower case spatial connections. We could now examine 
how spatial geodesics are deformed in the boost metric, 
but since these are not often used in practice we elect not 
to do so here. The standard expression, 

Ry = Ry - 2D[,C\]j + 

= Ry - 2D[,C\.], - , (33) 
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similarly relates the two curvatures. Note that in the 
three spatial dimensions of the spatial slice the Weyl ten¬ 
sor is identically zero, so we need only consider the Ricci 
tensors rather than the full spatial Riemann tensors. Pro¬ 
jecting upper case covariant derivatives into the lower 
case foliation immediately reveals the geometric mean¬ 
ing of the boost covariant derivative. Let Sa and Sab 
denote upper case spatial tensors, related in the stan¬ 
dard way cni) to their projections Sa and Sab into the 
lower case foliation. Then we have, 

_L = DiS + Wv, (CnS) , (34) 


These equations serve as a prototype when looking at the 
more complicated systems that follow. In particular the 
Lie derivative terms for tt differs from what one might 
naively expect. 

Gravitational field equations: We denote = BK® 

Vij. Moving now to write the field equations in terms of 
the boost metric we obtain, 

dtgij = -2 a Ky -f 2 a A(iUj) -k C(p-av)gij , (40) 

and after delicate use of the first hypersurface con¬ 
straint da. 


for a scalar S, and, 

+ Wv, (CnSj) - X\jSk , 

^DiSjk — -t- WVi (ZItV'^j/c) X ijSlk X ikSjl , 

(35) 


for the tensors, with, 

= [Kjvi - 2 + ViVjAi] . (36) 


The general pattern can be read off from these equa¬ 
tions. The boost covariant derivative is the part of the 
projected upper case derivative formed from lower case 
spatial derivatives and the boost vector. The remainder 
depends on the Lie derivative along and the boost 
extrinsic curvature Kij. We can relate the projected up¬ 
per case spatial Ricci tensor and the boost curvature by 
straightforward, albeit tedious, direct computation, 

R„- + DkX\j - D.X^k - X\iX\j 

+ X^uX^3k - Wvi T + Afi^^K) 

- W-^v'^ T 

- W-^v’^ T . (37) 


This result holds regardless of the gravitational field 
equations, but possible further manipulation is possible 
by the addition of the hypersurface constraints. The 
projected upper case covariant derivatives here can be 
replaced using (I55|) , which results in a slightly longer ex¬ 
pression in terms of IK^. 

Dual foliation formulation of the wave-equation: As 
a simple example we consider a 3 -I- 1 decomposition of 
the wave equation = 0 using the dual foliation. For 
the Lie-derivative of the boost metric we define. 


; (38) 


again with the convention that P = (g~^)*^ Pij. Intro¬ 
ducing the reduction variable /Catc/) = W-k. We then ob¬ 
tain. 


dtfi (y 7T — , 

dt^ = a (g-i)*^' + A, Bjfi] 

-f a [K -f Kyy -|- P -|- log(IF a)] tt -|- Z!(,g_|_Q,„)7r. 

(39) 


dt IKij = a Ry — W ^D(i [W A^)] -I- v'^a BKiOKj 

- W^a (DfcA/ + (2 - v^) Mi)v,vj 

- aC{yy-iy){w- 2W~'^D(^i\W‘^aKj)y\ 

- W^alAyKv.v, 

- 2 ex BK (^2 ^j)v T cr ([K -t- -t- P -t- Dy log(IF ck)^ 

o IK Pij — (X ^r?^K -j- 2 K^.„ -t- P -t- VF ^/JKy^ /JK(S:Vij 

- 2q;(K — A( 8 )i;)^’i (K — A®v)jfe — 2 Ky(^i 

- 4a(g"i)'='(K - A®-(;)fc(iPj), -f-Z:(^+„„)Ky . (41) 

Notice that we end up here with equations involving prin¬ 
cipal derivatives of the lower case shift but not the lower 
case lapse. That is, in equation dMl) first derivatives 
of /3* appear, but in equation (1411) no second derivatives 
of a appear, and instead we have derivatives of A^. This 
should be compared with the form of the equation (ED 
in which second derivatives of ex are present. Intuitively 
this happens because we are mixing the use of lower case 
spatial coordinates with upper case spatial objects. The 
Hamiltonian constraint can also be expressed in terms of 
these variables, giving, 

H = R -h (K -P K™)2 - Ky K*^' + 2Di (u* K) 

- 2 (g-i)*^ D,(K,„ + W-^'^ D[feWu,]) 

+ W-2 7 *^ 7 '=' D[iWvk] D[, Wu,] . (42) 

Likewise the momentum constraint becomes, 

nvi. = DJ{W-^Kk^) - D,(K + K™) -b z:,K„ 

+ W{g^^y'^[D,D[kWv,] - D,(IF-^' D[iWvk])] 

- 2(g-i)^S' D[iWv^-^ + R„ + P K„ , (43) 

where, up to additions of the hypersurface constraints, 
we define IH = iZ — 2XIy and DVIi = Mi. For readability in 
equations (USD and (I43p we write D^^AVvj-^ = D[i(IFfj]). 
No complications arise in including the stress-energy ten¬ 
sor for nonvacuum spacetimes. We note in passing that 
by taking g^ and K 7 as evolved variables the evolution 
equation for the metric also looks natural, since the pro¬ 
jected upper case acceleration appears as it would in a 
tetrad formulation when the timelike leg differs from the 
normal vector n“. With this choice of variables, we also 
see that the constraints can be written so that there is 
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no explicit appearance of the projected upper case accel¬ 
eration Ai. Furthermore, there is no explicit appearance 
of in equation dm)- Nevertheless we may be more 

interested in how the projected upper case extrinsic cur¬ 
vature evolves. This slightly more compact expression is 
trivially obtained by combining dm and m- The equa¬ 
tions are given in both forms in mathematica notebooks 
that accompany the paper (2^ . One expects to be able to 
formulate an initial data construction strategy naturally 
around the boost metric. A natural starting point for 
this would be to examine conformal flatness of the boost 
metric in boosted Schwarzschild. See [2^ for work along 
these lines. Observe that the notation in these last equa¬ 
tions is made slightly more cumbersome by continuing 
to insist on raising and lowering indices with the spatial 
metric jij , but we prefer to do so to avoid confusion with 
the surrounding calculations. 

The Generalized Jang equation: The Jang equa¬ 
tion is a quasilinear partial differential equation of 
minimal surface type, originally introduced as a tool for 
the proof of the positive energy theorem. Given initial 
data {'yij,Kij) for the initial value problem in General 
Relativity it reads, 

D^TD^T W D,D,T \ 

l + D>^TDkTj\ (l-hL»'=TiJfcr)i/2 j 

(44) 

for the unknown scalar field T. This equation was moti¬ 
vated by the characterization of slices of the Minkowski 
spacetime, in which there exists a scalar function T such 
that the second bracket of (HU) vanishes, and such that 
the boost metric, 

gy = lij + DiT DjT, 

is flat. With the present formalism it is clear that the 
natural curved space generalization to (HU) should be, 

(Ky - A(,u,)) = K - W-2 a„ = 0, (45) 

which, remarkably corresponds to the upper case foli¬ 
ation being maximal, because using the inverse boost 
metric to trace projected upper case quantities reveals 
the full trace of the original upper case tensor. One 
furthermore expects an analogous characterization of 
general asymptotically flat initial data sets to that of 
flat-space, in roughly the following terms: Gonsider 
data extracted from a spacetime {M,g), 

written in coordinates Xtt- = {T,X-), on some Gauchy 
slice, not necessarily a level set of T. An initial data 
set corresponds to the same data if and only 

if there exist vectors (u®, A^) satisfying the hypersurface 
constraint (HI, which we may write in the form, 

{D X Wv)i = {D\nA x A-^Wv)i , (46) 

and a projected Jacobian transformation ip-i with in¬ 
verse such that, 

^ij — Jij kF ViVj , 

K,, = - W-^D(iWv,) , (47) 



satisfy, 

'"'Vy = 3_m > 

- A(,u,)), (48) 

everywhere on the constant-t slice. The relationship be¬ 
tween the projected and true Jacobian transformation is 
that, 

gA, = + Wv^V,. 

(49) 

The transformation (f-i has the property that it maps N- 
spacelike contravariant tensor indices in X- coordinates 
to a:® coordinates and simultaneously projects the object 
into the lower case slice, and vice-versa for covariant n- 
spacelike indices in coordinates x®. The inverse property 
is easily verified by direct computation. The equivalent 
transformation <?®j in the opposite direction is defined 
in the obvious way. In the special case that the upper 
case coordinates are global inertial on Minkowski space¬ 
time this characterization reduces to that stated above 
motivating the Jang equation. We propose that this, 
rather than the conformal transformation, as is some¬ 
times claimed, constitutes the relation that should be 
used to build the natural equivalence class over phys¬ 
ically equivalent solutions to metric-based formulations 
of GR. 

Discussion of the dual foliation initial value problem: 
Given a boost metric gy, projected extrinsic curva¬ 
ture IKy, boost vector Vi, and acceleration Ai satisfy¬ 
ing the hypersurface, Hamiltonian and momentum con¬ 
straints we have a suitable set of initial data for vacuum 
GR. None of these quantities are invariant under changes 
of the upper case time coordinate T, but the form of 
the held equations is nevertheless invariant under this 
change. One way to view the resulting additional free¬ 
dom is that by breaking the correspondence between the 
spatial metric and the projection operator onto slices of 
the foliation, we gain the freedom to take the spatial met¬ 
ric of other foliations as the evolved variable. The sub¬ 
sequent projection into the foliation to obtain the boost 
metric is the most convenient way to deal with the vari¬ 
able in the 3-1-1 language, and Hts nicely with earlier 
work such as the Jang equation. It is natural to com¬ 
pare this reformulation with the freedom in the Maxwell 
equations, whose gauge can be altered without changing 
the coordinates on spacetime. The electric and magnetic 
fields are invariant under such changes, as they depend 
only on the Faraday tensor and the choice of coordinates. 
We have exactly the same status with the boost freedom; 
the electric and magnetic parts of the Weyl tensor are de¬ 
termined purely by the choice of coordinates, according 
to (H51) . and thus independent of the choice in the boost 
freedom. But the form of the field equations is invariant 
under changes to the boost. 

Working with the upper case spatial tensor basis: Al¬ 
lowing the boost freedom decouples, in the highest 
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derivatives, the lower case lapse from the evolved vari¬ 
ables. Therefore it is natural to ask whether such a de¬ 
coupling can also be obtained in the lower case shift by 
keeping all tensors in the X- coordinate basis. We there¬ 
fore now wish to drop the Jang-equation style use of two 
time coordinates with everything expressed in the coordi¬ 
nate basis a;* vectors, and instead use the X- basis tensor 
components, whilst computing derivatives in the co¬ 
ordinates. The time derivative of the projected Jacobian 
is, 

= £(/3-a«)¥>-i - a {AW)~^iphd^B-, 
dt{‘P~^yi = + 

( 50 ) 

Given the time derivative of an upper case spatial tensor 
in lower case coordinates we can now use ([HHI to compute 
the lower case time derivative in the upper case basis. 
Take for example Sat symmetric upper case spatial, again 
with projection Sij into the lower case foliation. Suppose 
we have, 

dt^ij — Gi Xij -j- B^Sij , (^1) 

then it follows that, 

dtS,^ = a + 2 a 

+ ( 52 ) 

Taking Sab as the upper case spatial metric, and looking 
at pni) we immediately see that indeed the lower case 
shift does become decoupled in the highest derivatives. 
It is sufficient to consider only this evolution equation 
because this is the only place where the shift is coupled 
in the principal part. The relationship between the upper 
case Christoffel symbol and that of the boost metric is, 

(<^-')\+ x\j 

+ - A-^Wv^,^,)di_B^ , ( 53 ) 

with X^ij defined as above (155)) . Note that dSSl) can be 
rewritten as, 

= x\, - A-^Wv^,^^-,)diB^, 

( 54 ) 

which in the setting with vanishing boost vector can be 
interpreted as the statement that the Levi-Civita connec¬ 
tion of the spatial metric is the gauge covariant deriva¬ 
tive associated with spatial diffeomorphisms. The full 
field equations are given in this mixed basis form in , 
together with a canonical Hamiltonian treatment. 

III. DOUBLE-NULL FORMULATION 

In this section we work in coordinates x^ = (t,a;*). 
Throughout latin indices a, b, c, d, e will be abstract 


as in the previous section, and likewise latin in¬ 
dices i, j, k, I, m,p stand for spatial components in x^ as 
before. We perform a 2 -|- I decomposition against r on 
the spatial slice, and take x^ = to be adapted 

coordinates. Thus upper case latin indices A, B, C, D 
stand for those in the level-sets of r. 


A. 2 -I- 1 + 1 Decomposition 

Motivation: We now turn our attention towards find¬ 
ing coordinates suitable for studying the collapse of grav¬ 
itational waves. It is known that sufficiently small per¬ 
turbations of the Minkowski spacetime are long-lived in 
the pure harmonic gauge [2 (tL l27j| , but this class of initial 
data presumably does not include every possible data set 
that eventually asymptote to the Minkowski spacetime; 
for sufficiently strong data, or indeed sufficiently strong 
pure gauge perturbations, coordinate singularities are ex¬ 
pected to form. Indeed there are examples of this phe¬ 
nomenon [ 2 ^ , and it may be that some of the difficulties 
in evolving strong Brill waves in [ 2 ^ were caused by the 
use of the pure harmonic slicing. Therefore we look else¬ 
where. Empirically the generalized harmonic gauges [s^ 
have been found very robust in both binary blackhole 
and collapse scenarios. However, from the mathematical 
point of view, the strongest results concerning collapse to 
a blackhole employ a double-null foliation Q . It is known 
that a particular type of initial data will form an appar¬ 
ent horizon before any coordinate singularity forms. It 
is not clear that close to the critical threshold of black¬ 
hole formation these coordinates are well-behaved. It is 
also rather doubtful that the double-null foliation will be 
useful in the strong-held region in binary-blackhole space- 
times. But since these coordinates naturally conform to 
the causal structure of the spacetime, and there is like¬ 
wise no guarantee of nice behavior for any other coordi¬ 
nate system, they seem to be in the best shape for con¬ 
sideration. In particular, in both 3 -|- I-spherical symme¬ 
try [3l| and a 2 -|- 1 dimensional setting |3^ such coordi¬ 
nates have been effectively used in studies of critical col¬ 
lapse, neatly sidestepping the need for mesh-rehnement. 
We thus look at related gauges suitable for the initial 
value problem. 

The 2 -|- 1 decomposition: In a double-null foliation 
there are two crucial coordinates, optical functions, 
whose level sets are incoming and outgoing null surfaces. 
In the 3-1-1 setting we have however only singled out 
the time-coordinate for special treatment. A given pair 
of these null hypersurfaces intersect in a spacelike two- 
sphere. Given a spatial slice of constant t, complete with 
spatial metric and extrinsic curvature Kij let us de- 
hne a new coordinate r, which we will use to perform 
a 2 -I- 1 split. The idea is that the level sets of r should 
become the spheres on which the null surfaces intersect. 
Obviously the calculations that follow are essentially the 
same as those of the standard 3-1-1 split, as described in 
detail in [T|. The coordinate r defines a unit normal s* 


to a surface of constant r according to, 

L-2 = j^i{D,r){Djr ), s* = . (55) 

We will call L the length scalar. The normal vec¬ 
tor s® naturally defines the induced metric in the two- 
dimensional level set, 

qij = lij - s^SJ . (56) 

Likewise we have the extrinsic curvature, 

Xzj = q'^^DkSj , (57) 

so the first derivatives of the spatial normal vector are in 
total, 

DiSj = Xij - Sipj InL . (58) 

We denote the covariant derivative compatible with the 
induced metric qij by pi^ and likewise use for the par¬ 
tial derivative projected into the surface. Note the rel¬ 
ative change in sign in the definition of this extrinsic 
curvature and that of the slice Kij. Let the vector r® be 
tangent to lines of constant 0"^, the two spatial coordi¬ 
nates in the level set. We have 

r® = Ls® -f 6®, (59) 

with Vsi = 0. We call the two-dimensional vector 5® 
the slip vector. Note the relation r^'Dir = 1. With this 
notation we can express the spatial metric as, 

df = -f qAB {de^ + h^dr) (d6»® -f b^dr). (60) 

When performing the 3-1-1 decomposition one finds that 
the lapse scalar and shift vector are freely specifiable, 
which is of course not the case with the analogous length 
scalar and slip vector. The four-dimensional metric is, 

ds2 = -a^dP + L‘^{dr + L-^^^dtf 

-f qAB{de^ + h^dr + /3^dt)(d6»® -h h^dr + P^dt). 

(61) 

The coordinate light speeds in the increasing and decreas¬ 
ing r directions are, 

ci = (-/3* ± a) L-i, (62) 

whilst in the transverse directions we have, 

ci = -l3^Tb^aL-^. (63) 

Although they may not necessarily be associated with 
spherical-polar coordinates we will call these transverse 
directions ‘angular’. An obvious choice is a = L, 
and /?'’ = 0, under which the coordinate light-speeds 
are c!j_ = ±1. With this choice the combinations u = t — r 
and V = t + r are the optical outgoing and incoming null 
coordinates alluded to earlier, and the coordinates are 
naturally adapted to the causal structure of the space¬ 
time in the null ± s'® directions. 


The extrinsic curvatures: We immediately split the 
extrinsic curvature xab into a trace and tracefree part, 

XAB = XAB + ^AB X ■ (64) 

We will similarly use the notation q for the determinant 
of the two-metric. Finally we decompose the extrinsic 
curvature of the spacelike surface as embedded in the 
spacetime by 

ATjj — SiSjKss T 4” ‘^^{iq j)^A d~ Kab ; (65) 

where in accordance with the previous notation Kab 
stands for the projected, tracefree part. We use in¬ 
dices s to denote contraction with s'®, and qq for a 
trace taken with the two-dimensional metric qAB- We 
write Kab = q^Aq^BKij for the projected extrinsic cur¬ 
vature, and occasionally still use K = Kgs + Kqq as a 
shorthand for the trace of the extrinsic curvature. 

The Christojfel symbol: The spatial Christoffel sym¬ 
bol is readily decomposed as, 

r^ss=Cs{\nL), r%,g®A = ^A(lnL), 

(?®A q^B = —XAB , q'’A q’B qk^ =J^AB, 

T\gq^k = -^(InL) + ^ {Crb^ - b^^Bh^) , 
T’‘^s q^B qkA = XAB + \qAC^Bb^ , (66) 

with ,r denoting the Christoffel symbols of the two- 
metric qAB- The spatial contracted Christoffel sym¬ 
bols r® = are therefore, 

T^ = Cg{\nL)-x, 

r®g/ = + {CM - b^^Bb^) - lfA{\nL). (67) 

Curvature: As in (16511 the spatial Ricci curvature Rij 
can be decomposed according to, 

Rss = -jjp'^pAL - CsX - X^BX^A , 

Rqq = -tP'^V>aL - CsX + ^- 
RsA = pBX^A - Pax ; 

Rab = -CsXab - ^PaPb^L + 2x^aXbc , (68) 

by the classical Gauss-Codazzi-Mainardi equations. 

Vacuum field equations: The 2 -|- 1 -|- 1 decomposed 
Einstein equations are given first by the constraints, 

H = fi- ipAp^h - 2CsX - XabX^^ - lx" 

- iK^qq - 2KqqKss + 2 KaK^ + KabK^^ , 

Ms = -CsKqq + PaK^ + 2K^pA{\nL) + xKse 

- \xKqq+X^BK^A, 

Ma = CsKa + P^KaB — ipAKqq — PaKss + XKa 

- KsspA{\kvL) + \KqqpA{\nL) + kABp^{\kiL). 

(69) 

Since null geodesic expansions in the n‘®±s‘® directions are 
given by x± AT^q, we can view the H and Mg constraints 
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as dictating how the expansions vary over the slice. Next 
we have evolution equations for the metric, 

^((InL) = -aKss + + 0^(3“ , 

dth^ = -2aLK^ + L^Ijf{L-^l3p + £^/3^ + ^^ 36 ^ , 
dtqAB = —2aKAB + 2(3 ^xab + If-pqAB , (70) 

and for the decomposed extrinsic curvature, 

dtK,, = -CsCsa + a[Rss + 2KaK^ + 

- IfP{\AL)pAa - 2LK,ML-^pp 

- 2LK^pA{L-^n + P^DsKss + P^PaKss , 

dtKqq = -JjPpACc + + Kl^ + 2KaK^ + 

- xCsOc + 2LK^pA{L-^n + 

+ P^pAKqq , 

dtKA = —pA^^sOi + 0:[RsA + KqqKA] + ApBOi 

- LKADs{L-^n - LK^ApB{L-^n 

+ (f-P^A + C-sKa ) 

dtKAB = -pApVa - XABC.a + 2LK(^ApBf {R-^PP 

+ o:[Rab — “^K^aKbc + (KaKb) + KssKab] 
+ P^CsKab + 4-pKab ■ (71) 

Here we have defined the Lie-derivative in the level-set 
of r, in the obvious way, and where it is understood that 
the vector argument must be projected with qAB, which 
allows us to write, for example, 

4pKab = S ( qPP^PjKu + ^k^^^Pr i<iPP')) ■ 

(72) 

Finally the equation of motion for the extrinsic curva¬ 
ture Xij can be computed from the relation 2xij = Csqij- 
The result is, 

dtXAB = —J^sicxKAB) + aJ^KqAB + ‘2 K(^aPb)<A 

+ OiKssXAB + ‘^OiK (^j^Pb){}ai L) -|- {j-pApBL)P^ 

~ PaPbP^ + P^^sXAB + ^PXAB ■ (73) 

If we want to treat r as a radial coordinate, and the 
remaining 9^ as angular coordinates we obtain regularity 
conditions at r = 0. These conditions will be discussed 
elsewhere. 


B. Double-null formulation 

We now look at the field equations imposing the 
double-null gauge explicitly. The aim here is first, to 
present the simplified form of the field equations in this 
gauge, and second, to examine whether or not hyperbol- 
icity of the full system can be obtained. 


Pure gauge analysis: Let us examine the behavior 
of infinitesimal perturbations to coordinates satisfying 
the optical conditions a = L and /3® = 0 above, ad¬ 
ditionally taking = b^. This sign is chosen as¬ 
suming that the gravitational wave is traveling mostly 
in the minus r direction consistent with earlier work. 
This is Christodolou’s gauge choice in rewritten in a 
time-space rather than a double-null form. The expres¬ 
sions |10l | for the time development of the perturbations 
to the time and space coordinates are, 

dte = U- + P^diO , 

dtip^ = V^ + aDP - 6DP + Cp'tp^. (74) 

Then we have, 

U = A[a] = A[L] = DirDjr)-^/'^] 

= LCsP-OLKss+VPaL, 

= s,A[P ], 

= A[P^] = A[b^] 

= -2eLK^ + l2^(L-V«) + Cri:^ + 4.^b^ , 

(75) 

and obtain the pure gauge subsystem, 

dt[L-p] = LZ4,[L-Vs] + b^pA[L-p] , 
at[L-Vs] = LDs[L-p] + b^pA[L-^i^s \, 
dt{q-ip)^ = dr{q-+ L'^P^[L~'^{9 PiIjs)] 

-2LK^{e + ijs)- (76) 

This first order PDF system is only weakly hyperbolic. 
The arguments presented in 0, building on those 
of [§, [13 , can be used to show that no strongly hyper¬ 
bolic formulation can be built with this gauge condition, 
at least if the formulation is constructed under the stan¬ 
dard free-evolution approach. Therefore the double-null 
gauge can not be directly used in numerical relativity in 
the standard way, but requires a more subtle approach, 
or some modification. Note that the problem here comes 
from the choice P^ = 6"^, and there are simple modifica¬ 
tions under which strong hyperbolicity of the pure gauge 
subsystem can be obtained. Nevertheless, building a for¬ 
mulation of GR which is at least strongly hyperbolic with 
one of these good, modified, conditions will be more in¬ 
volved because here the s* direction is singled out for 
special treatment. Therefore in the following sections we 
instead look for a simpler approach employing the dual 
foliation formalism. 

Fixing the gauge: Despite the shortcomings un¬ 
earthed by the pure gauge analysis, for completeness we 
present the full held equations with in the double-null 
form. As above we choose a = L, /3® = 0 and P^ = b^. 
This choice has no effect on the constraints (|551) . but the 
evolution equations become, 

dt{\nL) = -LKss + b^pA{\aL ), 
dtb^ = -2L^K^ + Crb^ , 
dtqAB = —2aKAB + 4bqAB , (77) 
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for the metric components, and 

dtK,, = -CsLsL + L[R,s + 2 KaK^ + + KI] 

-lp^{\nL)]pAL + b'^]pAKss, 

dtKqq = -JjfjpAL + L{Rqq + + 2KaK^ + 

- xCsL + b^pAKqq , 

QiKa = ~1PaC-sL + L[Rsa + KqqKA] + A JpsL 

+ iJ-pKA , 

dtKAB = —IpApB^L — xabCsL + L[Rab — 2K^ aKbc 
+ {KaKb) + K.sKab] + tlbKAB , (78) 

for the extrinsic curvature. Finally we have, 

dtXAB = —C,s{LKab) + LI^KqAB + 2K(^aPb)L 

+ LKssXab + ‘2LK(^a^b) (111 7 j ) + J^bXAB ■ 

(79) 

Taking linear combinations of these variables one can 
rewrite so that all of the equations take the form of 
‘transport equations’ in the n“ ± s“ directions, but with 
transverse derivatives appearing as sources. Particularly 
relevant are the combinations x ^ Kqq, as they reveal 
the Raychaudhuri equations. Up to this trivial change of 
variables and the split into time-space derivatives, rather 
than the double-null choice, this is the same system pre¬ 
sented in Ch. 3 of 0, where it was also noted that this 
system is not hyperbolic. 

IV. COORDINATE SWITCHED FIRST ORDER 
GENERALIZED HARMONIC GAUGE 

A. Double-null Jacobians 

Time and Radial coordinates: Let us now abandon 
the idea of evolving in double-null coordinates directly, 
and instead examine how the spacetime could be con¬ 
structed in these coordinates a posteriori, having con¬ 
structed the spacetime locally in the harmonic gauge, for 
example. This is similar to the strategy employed in Ja. 
It has also been used in numerically, in for example [34| . 
Let us work in lower case coordinates x^. We would like 
the upper case coordinates to satisfy the double-null con¬ 
ditions. As elsewhere, 

= -A V“T, S"* = L . (80) 

with AL = (T, R, 0—). We choose A = L to impose the 
double-null gauge, regardless of the angular coordinates. 
Under this condition we define ingoing and outgoing null 
vectors L“, K°‘, 

-N‘^ - S‘^ = K‘^ = LK^. 

(81) 

The renormalized vectors and 741“ generate null 
geodesics in the and 741“ directions. It is natural to 


define two lower case spatial vectors and from the 
Jacobian according to, 

vf = E±sf = ^,T L-^Wv,. (82) 

The vectors have unit magnitude. Indices s± stand 
for contraction with these vectors. The scalars E± are 
the energy of the ingoing and outgoing congruences as 
measured by the Eulerian observers n“. In terms of the 
Jacobian they are, 

E± = L-^W T ■ (83) 

The null geodesic vectors are then, 

7tl“ = -E_ (n“ -b s“ ) , ii“ = -E+ (n“ - s“ ) . (84) 

Straightforward computation then reveals evolution 
equations, 

dt In E± = C(^i3±as±) \nE±+a {Ks±s± ± In a) , 
dtvf = ±a ± E± D^a + Cpvf . (85) 

The hypersurface constraints were used freely to arrive at 
this result. Interestingly the term appearing as a source 
in the first equation is essentially a characteristic variable 
of the (first order in time, second order in space) GHG 
formulation. Notice furthermore that, after adjusting the 
normalization of these can be compared with the re¬ 
sults of and are of course compatible. Note also 
that the equations for E± follow from = E\. 

The equations (1851) form a symmetric hyperbolic system 
in {E±^v^) which is equivalent to a system in the scalar 
components L~^W and tt— and the vector parts (f>E^ 
and L~^Wvi of the Jacobian. 

Angular coordinates: To complete the equations of 
motion for the Jacobian we require a choice for (p^i 
and TT—. We have already seen that from the pure gauge 
point of view the choice ‘/3"^ = b"^' is problematic. To 
understand how this issue appears in the Jacobian for¬ 
mulation, let us assume that the component tt— takes the 
form, 

TT— = , (86) 

for some known lower case spatial vector m®. This cor¬ 
responds to determining the angular coordinates by Lie¬ 
dragging so that (n“ -I- m“)Va0— = 0. This family in¬ 
cludes = ±6^’ by choosing m“ = ±s^. Plugging the 
relation (1551) into the Jacobian equation of motion (1151) 
and using the hypersurface constraints gives the compact 
expression, 

dt(jAi= C(^p_am)4^i- (87) 

If the vector m“ is given a priori this equation is hy¬ 
perbolic. On the other hand if we wish to make m“ 
dependent on the other components of the Jacobian, 
the derivative of m“ in the Lie-derivative is problem¬ 
atic, as it is a one-way coupling, in the sense that the 
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equations of motion for the angular coordinates explic¬ 
itly depend upon the (T, R) coordinates in the princi¬ 
pal part, but not vice-versa. This type of coupling is 
dangerous because it can leave non-trivial Jordan blocks 
in the principal symbol if the speeds associated with 
the T, R and 0— blocks clash. Fortunately this dis¬ 
cussion suggests two alternatives. The first is to con¬ 
struct the angular coordinates using just vectors associ¬ 
ated with lower case coordinates. For example we could 
choose n“Va0— = 0, or even t“Va0— = 0, in which case 
the angular coordinates could correspond to those built 
directly from the lower case coordinates. One possibil¬ 
ity would be to use a reference metric to build a first 
order GHG formulation directly in spherical-polar coor¬ 
dinates. Then the upper case angular coordinates could 
be exactly those of the lower case system. The second 
option is to make sure that the speeds of the two sub¬ 
systems do not coincide, or that if they do the princi¬ 
pal symbol is nevertheless diagonalizable. For this we 
might try iV“Vo0— = W (n“ -I- w®)Va©— = 0, or in other 
words m“ = z;“, which is also equivalent to B— = 0 when 
working in upper case coordinates. Since, 

= iE++E_)-^{vl-vl), (88) 

we arrive at, 

= —av^Dj4>^i + ^aLW~^(tAjDi{v^j^ — v^_) 

- \aLW-^ t^-D,{E+ + E_) + TT^Aa + Cp(lAi , 

(89) 

after expanding the Lie-derivative. 

Hyperbolicity of the Double-Null Jacobians: The 
double-null Jacobian system with to* given a priori is 
trivially symmetric hyperbolic. Choosing instead to* = 
u* as in (ESI), the principal symbol for the subsystem with 
the components (£1±,(/e^s) is, 

P* = aA'*-f/3n, (90) 

with, 

0 0 0 \ 

0 0 

0 ±4 0 

(91) 

with an index ‘s’ denoting contraction with an arbitrary 
lower case unit spatial vector s* and Pab = lab — SaSb 
as elsewhere. The remaining block of the full principal 
symbol, associated with is decoupled, and has no 
affect on the following discussion. The eigenvalues of the 
principal symbol are ,5® ± asj- and /3® — av^. For strong 
hyperbolicity we need that the symbol is diagonalizable 
for every s*. Choose for example s* such that, 

{E - s*_)s, = 0, (92) 


then the ‘sf_’ and ‘u®’ eigenvalues coincide, and the prin¬ 
cipal symbol is missing eigenvectors. Therefore the sys¬ 
tem is again only weakly hyperbolic. This type of degen¬ 
eracy occurs if TO* is taken to be any vector constructed 
from v\., thus we must abandon the second alternative 
suggested by the above discussion. This leaves the option 
to fix TO* a priori, or to choose a gauge of a different form 
for the angular coordinates, such as the generalized har¬ 
monic option da- At least for the choice = b^' this 
result was expected, because the pure gauge subsystem 
we examined before is closely related to the system for 
the Jacobians. 


B. First order generalized harmonic gauge 

There is much work about formulations of general rel¬ 
ativity in first order form. For a review, see [s^. Of 
particular interest in recent years has been the first or¬ 
der reduction of the GHG formulation [H, [13] • To a 
large extent this interest was driven by use inside the 
SpEC numerical relativity code. Here we summarize just 
the relevant features. The first order GHG system is a 
quasilinear, first order symmetric hyperbolic system of 
the form, 

9tU = A^c^pU -I- S . (93) 

The formulation has a set of constraints compatible with 
the evolution equations. The key feature of the system 
is that the principal part takes the form of a first or¬ 
der reduction of the wave equation. This is attained by 
carefully coupling the coordinate choice Dx" = to 
the held equations [l^ and then reducing to Hrst order. 
Gravitational radiation controlling, constraint preserving 
boundary conditions for the system have been studied 
and implemented [l^, 113, • The evolution system is 

now used frequently in the evolution of compact binary 
spacetimes. Of crucial importance in the present work 
is that in the hrst order reduction, all hrst derivatives 
of the metric can be expressed in terms of evolved vari¬ 
ables without taking any derivatives. So, for example if 
we were to evolve the double-null Jacobians (IRH|) along¬ 
side the GHG variables, we can formulate the system so 
that no coupling occurs through derivatives. In other 
words the Jacobians are minimally coupled to the GHG 
formulation. 


C. Coordinate switch and hyperbolicity 

Coordinate switch: Suppose, without loss of general¬ 
ity, that we are given a hrst order quasilinear system of 
the form, 

dTU = {AAP + BPl)dpU + AS. ( 94 ) 

Now consider the effect of a change of independent co¬ 
ordinates on the whole system, so that rather than us¬ 
ing the upper case Xh- coordinates we employ the lower 
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case ones. In the new coordinates the system of course 
retains the same functional form, but now with, 

{l+A^)dtu = aW-^ (1 + A^)np)5pu 

+ aW-^S. (95) 

Recall here the various components of the inverse Jaco¬ 
bian defined by m and (l4^ . and the shorthand A— = 
A-Vi- It is easily confirmed that any closed constraint 
subsystem remains closed. 

Symmetric hyperbolicity: The original system is as¬ 
sumed to be symmetric hyperbolic, that is, there exists 
a symmetric positive definite symmetrizer H such that 

H {AAP + BE-l) , 

is symmetric for each p. In the present context, this 
system can be taken as a first order generalized har¬ 
monic formulation coupled, minimally, to the double-null 
Jacobians (ESI). Under a smallness assumption on Vi 
symmetric hyperbolicity is unaffected by the change of 
coordinates. Taking the system in the form (IM|) . the 
symmetrizer is unaffected by the change of coordinates. 
If we insist on multiplying on the left by the inverse 
of (l -f A—), then we pick up exactly a factor of this 
matrix on the right of the modified symmetrizer. But 
to obtain energy estimates for metric components using 
norms formed from the lower case coordinate basis com¬ 
ponents of tensors more effort may be necessary. 

Strong hyperbolicity: The principal symbol of the sys¬ 
tem in Alt is, 

P|=AA^+R^1, (96) 

where on the right hand side superscript S_ denotes con¬ 
traction with an arbitrary unit upper case spatial covec¬ 
tor Si- Multiplying (IMl) by the inverse of (l + A—), the 
principal symbol in can be read off, 

V% = a W-^ ((1 + A^) A^ - nn) . (97) 

On the right hand side superscript s denotes contraction 
with an arbitrary unit lower case spatial covector Si, and 
underlined superscript s denotes the same contraction, 
but pushed through the transformation ((p“^)*i. Strong 
hyperbolicity is the requirement that for each Si there ex¬ 
ists a symmetrizer H®, uniformly symmetric positive def¬ 
inite in Si, such that the product of the symmetrizer with 
the principal symbol is symmetric. For quasilinear prob¬ 
lems strong hyperbolicity, together with some smooth¬ 
ness conditions on the symmetrizer are sufficient to guar¬ 
antee local well-posedness of the initial value problem. 
These smoothness conditions are sometimes included in 
the definition of strong hyperbolicity. The existence of 
a symmetrizer for fixed Si is equivalent to the require¬ 
ment that the principal symbol has a complete set of 
eigenvectors with real eigenvalues. Assuming that the 
system is strongly hyperbolic in the upper case coordi¬ 
nates Xtt this necessary condition can be seen to hold in 


lower case coordinates x^ as follows. This discussion is 
adapted from [^ . Fix s^. The characteristic polynomial 
of has real eigenvalues, and thus we have a hyper¬ 
bolic polynomial with respect to N^. It follows that if 
the boost velocity is sufficiently small then P| also has 
real eigenvalues [^. Estimates of the range over which 
this condition is satisfied can be given. In the presence 
of a simple block structure of A-, the eigenvalues will 
simply be multiplied by those of (1 -I- A—). Take one of 
these eigenvalues A. Strong hyperbolicity in AE implies 
the existence of symmetric positive definite H with, 

HL = H(A^-(l-f-Ai^)(A-fn®)) , (98) 

symmetric. Suppose that u is a non-vanishing vector 
in the nullspace of L, and thus either an eigenvector or 
generalized eigenvector of P® . Suppose that it is a gen¬ 
eralized eigenvector, so that u = Lv for some v. Then, 

u^Hu = u^HLv = v^HLu = 0. (99) 

The second equality holds by symmetry of H L and the 
third because Lu = 0. Since H is symmetric positive 
definite we then have that u = 0. Therefore if u is in the 
nullspace of L it is a true eigenvector of P®. Thus P® 
has a complete set of eigenvectors. Details concerning 
the remaining uniformity condition in Si can be found 
in [40l |. 

Discussion: The advantage of the dual foliation is 
now clear. In section IIII Bl we were unable to find a 
hyperbolic formulation using the double-null coordinates 
directly. Using the dual foliation however the construc¬ 
tion of such a formulation is essentially trivial. Since 
the whole system is symmetric hyperbolic we do not 
lose regularity when mapping between the two coordi¬ 
nate systems. Regularity will however need to be looked 
at carefully if we are to use the double-null radial coor¬ 
dinate all the way to the origin, but this issue is only 
that of using spherical polar coordinates, and not related 
to the dual foliation strategy. With a little more book¬ 
keeping we can instead work from a standard first order 
in time, second order in space [1^,143 formulation of GR 
and avoid the first order reduction, but postpone any 
presentation thereof. One expects that this formulation 
could be treated according to Q to give an economical 
demonstration of a ‘neighborhood theorem’ for the GR 
characteristic initial value problem. Physically what has 
happened is that the coordinate and gauge degrees of 
freedom GR were decoupled as much as possible. Note 
that in earlier work a symmetric hyperbolic formulation 
using a double-null gauge was constructed using frame 
variables coupled to the Bianchi equations [d^. In con¬ 
trast here the aim was to arrive at such a formulation 
with as few changes as possible to standard formulations 
used in numerical work. From this practical point of view 
there is always the danger that the matrix (l -|- A—) be¬ 
comes singular. We can mitigate against this by simply 
evolving for a short but finite coordinate time, and then 
resetting the Jacobian to the identity by transforming 
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all the fields into the upper case tensor basis. For sufii- 
ciently regular data hyperbolicity guarantees short time 
existence so that this procedure can be performed iter¬ 
atively. Coupled with strong theorems on continuation 
of solutions, this strategy could perhaps even be used 
to guarantee that numerical calculations progress suc¬ 
cessfully into extreme regions of spacetime, although the 
double-null gauge might have to be replaced with some 
other suitable choice. For the numerical relativist, proba¬ 
bly the simplest summary of the double-null Jacobian re¬ 
sult is that it is t he g eneralization of the dual coordinate 
frames approach |46j employed in the SpEC code [l^ to 
the situation in which two differing time coordinates are 
considered, and where the Jacobians satisfy dynamical 
equations rather than arising as derivatives of algebraic 
relationships between the coordinates. 

V. CONCLUSION 

Making a dual foliation approach, we have shown that 
it is possible to effectively decouple the choice of coor¬ 
dinates from local well-posedness of the field equations 
of general relativity. This was done by evolving a first 
order reduction of the generalized harmonic formulation 
alongside Jacobians mapping from the desired coordinate 
system to the generalized harmonic one. The important 
example of the double-null gauge was considered. But 
in fact the set of coordinate choices resulting in equa¬ 
tions of motion for Jacobians that are minimally coupled 
to the remaining field equations is extremely large, and 
local well-posedness inside this class follows from well- 
posedness of the coordinate choice alone. It is thus ex¬ 
pected that with due care, this observation will allow us 
to evolve the generalized harmonic formulation using say, 
the Maximal slicing, and quasi-isotropic spatial coordi¬ 
nates, for example. We hope that this would allow for 


direct comparison of the results of (4^ in a modern nu¬ 
merical relativity code with minimal changes. See 
for recent work on the problem. It is important to real¬ 
ize that the role of the generalized harmonic formulation 
in the construction is completely auxiliary. In fact any 
well-behaved formulation of GR is amenable to the same 
trick, perhaps after a suitable reduction to first order, if 
we wish to avoid complicated book-keeping. The general¬ 
ized harmonic formulation is simply the most convenient 
example. 

For numerical applications it will be necessary to 
translate the constraint preserving boundary condi¬ 
tions [ 5 ^ of the harmonic system into the new set of 
variables and coordinates, but no major difficulty is 
expected in doing so. The obvious next step is to sys¬ 
tematically implement and test the approach, preferably 
in a simple context. Afterwards we expect to use the 
double-null coordinates to study the critical collapse of 
gravitational waves using the bamps code [itI - Ii^ . 
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